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Using a recently introduced mapping between a scalar elastic network tethered at its bound- 
aries and a diffusion problem with permanent traps, we study various vibrational properties of 
progressively tethered disordered fractals. Different scaling forms are proposed for different types 
of boundary tethering and numerical investigation of the scaling is performed by the approximate 
diagonalization of the corresponding large, sparse transition probability matrices. Rather different 
localization behaviors are found for the leading modes depending on the types of tethering. 



I. INTRODUCTION 

In this work we consider the vibrational spectrum associated with the scalar elasticity of fractals when their 
boundaries are progressively tethered to immobile, external anchors. By scalar elasticity, we mean the vibrational 
problem where the local displacement variable Ui{t) at site i is a scalar. We study the various signatures of the 
vibrational density of states and the localization properties of the normal modes under these conditions. 

The spatial scale invariance of fractals |l|] and the absence of translational invariance due to the boundaries conspire 
to influence their vibrational spectra in fundamental ways. The now well-studied fracton H] spectrum, with the 
characteristic, fractional power law for the low-energy density of states, arises for the bulk fractals with free boundaries, 
while so-called fractino [0 spectrum arises for objects which do not have fractal bulk but have fractal boundaries 
that are clamped. Here, we generalize these further to study bulk fractals whose fractal boundaries are progressively 
clamped or tethered to immobile anchors. The resulting vibrational spectrum is complex yet shows some simple 
scaling features common to many phenomena with long range order. 

The problem of vibrations in a disordered system is in itself of physical interest. Granular systems Ej, such as sand 
or snow piles or crops stored in silos, are among the systems of potential application. For example, sound propagation 
in a granular system had been modeled by Leibig pi using a scalar elastic network with a random distribution of spring 
constants and analyzed in terms of the normal modes. Boundary conditions such as tethering add much complexity 
to the problem which may be relevant to some physical and engineering applications. For example, random mixtures 
of soft and hard materials such as solid granules embedded in a polymer gel or colloidal particles and aggregates 
filling the pores or cracks in rocks may constitute such an application. Another potential interest might be due to the 
recently discovered anomalous behavior [pj of the superfluid transition when aerogels H are embedded in liquid "^He 
or '^He, where the aerogels may impose a boundary condition over a fractal boundary of the liquid helium. 

Most previous works on the effect of clamped boundaries on the vibrational density of states of inhomogeneous 
systems dealt with deterministic inhomogeneities (such as in fractal drums 1^). In this case the bulk is nonfractal 
(Euclidean) and the boundary is an ordered fractal, and the effect is a correction term to the leading behavior, 
the latter remaining the same as for a homogeneous system as one approaches the asymptotic limit of large system 
sizes. We have a substantially different system where both the bulk and the tethered boundaries are statistical 
fractals, sometimes of different fractal dimensionality, where even the leading behavior is expected to be special to 
inhomogeneous systems. 

The prototype of such a system is a critical percolation cluster Q . We have used site percolation clusters created 
near the critical percolation threshold on square and simple cubic lattices. We focus our attention on the low-energy 
regime where we expect the clamping effect of the fractal boundaries to be pronounced. Without tethering, these 
modes tend to have relatively large spatial extent, and thus are more sensitive to boundary tethering. In particular, 
we study the lowest energy mode (ei) and the peak (e^) that appears in the low-energy region of the vibrational 
density of states. The lowest energy mode is computationally much easier to obtain than the entire density of states 
and yet captures some of its essential features. 

To study these systems, we exploit the mapping between scalar elasticity and diffusion as extended |p[ to the case 
of mapping between the vibration of tethered objects and diffusion with permanent traps. In this approach, we define 
a transition probability matrix W (an S x S matrix for a cluster of S sites) where Wij is the hopping probability per 
step pij from site j to site i in the corresponding diffusion problem. We set pij = 1/z for all pairs of nontethered sites 
i,j where z is the lattice coordination number. The diagonal elements Wu = if site i is tethered (representing the 
complete leakage of diffusion field or /w/Z tethering) while, if i is not tethered, Wu = l — Ui/z (representing conservation 
of diffusion field) where rii is the number of available neighbors for site i. Then the time evolution of the diffusion 
field Pi{t) is given by 

P,{t + 1) = Y.[p^,P,{t) + {S,, - p,,)P{t)] = Y^W^.P.it) , (1) 

J 3 

where ^ includes the diagonal terms. 

The cigcnmodcs of W are the normal modes of the tethered vibration problem and the eigenvalues A are related 
to the classical vibration frequency uj and energy e = o;^ by 

A = 1 - cj2 « exp(-tj2) (2) 

where the last approximation corresponds to the long time limit. We use this formulation in part to take advantage of 
the fact that the low-energy region of the spectrum corresponds to the region of the maximum eigenvalues of W [lC| ] 
which affords a much easier access numerically. The numerical technique used is based on the work of Saad |T^] which 
implemented Arnoldi's method. Preliminary results of this work were reported earlier in pj which introduced the 



mapping as well as presented preliminary numerical results only for the all boundary tethered case in two dimensions. 
In the current work, we have improved data in both two and three dimensions as well as for hull tethering (see below). 

II. SCALING OF LOWEST ENERGY MODE FOR HULL VERSUS ALL BOUNDARY TETHERING 

In a cluster with tethered boundaries the lowest energy mode ei crucially depends on how many and which sites are 
tethered. Clearly, for an individual untethered mode, tethering of the sites where the mode has a large component 
greatly affects the mode. Statistically, an important classification of the boundary sites is between the internal and 
external boundaries (or hull) [|2|. For a fractal like the critical percolation cluster, they are both fractal objects with 
similar (but not identical) fractal dimensions; in two dimensions, the hull has a fractal dimension d'l = 1.75 while 
the internal boundaries have a fractal dimension d^r — 91/48 (same as the bulk fractal itself and also as the fractal 
dimension of all boundary dt). However, the behavior of ei (and the density of states) is unexpectedly different 
depending on if only (some of) the hull sites are tethered or (some of) both types of boundary sites (i.e., the hull and 
internal boundary) are tethered, as functions of the number N of the tethered sites and the size S of the cluster, the 
two main parameters in these situations. 

In order to understand the interplay of the different length scales we write the scaling variable as the ratio of the 
relevant length scales for the two cases. The common length scale for both cases is the average distance I between the 
tethered sites where 

l^iS/N^/'^f. (3) 

We also extend the meaning of A^ by TV = fNr where / is the fraction of the Nt boundary sites which are tethered: 
Nt = Nh (number of hull sites) for hull tethering, Nt = Na (number of all boundary sites) for all boundary tethering. 
In particular, this means N can be smaller than one if / is sufficiently small, which indicates many random samples 
in the ensemble will have no tethers. 

The second relevant length scale L for hull tethering for the critical percolation cluster is the the average diameter 
of the cluster itself, 

Lc - S^/'^' (4) 

since the hull forms a connected set of external boundaries to whose interior the normal modes are confined. In 
particular, the spatial extent of the lowest energy mode is dictated by that of the interior sites at the scale of the 
cluster diameter L^)- See the illustration of these lengths in Fig.l. 

On the other hand, in the case of all boundary tethering, the second relevant length Lf, is the average diameter of a 
blob in the cluster. As sites of both hull and internal boundary are tethered in all boundary tethering, a normal mode 
is confined to sites which are internal to both boundaries, i.e., in the blobs. The blobs in this context may be the 
areas of the cluster with essentially no internal holes (empty sites in the interior), which are connected together by 
links of less connectivity to form the overall cluster. See the illustration of these lengths in Fig.l. Though we do not 
have a precise characterization of a blob, it is clear that this sense of a blob is different from those discussed earlier, 
e.g., in the sense of a merely multiply connected component |l3[ . 

Since one would expect these new blobs to go critical at the ordinary critical percolation threshold pc, its diameter 
Lf, may be expected to follow a power law at Pc, 

Lb ^ Sy (5) 

for an appropriate universal exponent y in the scaling regime of large S. We note that both a precise identification of 
such blobs and the relation between Lb and S are still open questions |14[ |. 

Let us consider scaling for the hull tethering case first. In the limit of no tethering, the problem is reduced to 
that of so-called ants Q where the lowest energy mode is the stationary mode with e = for any S. On the other 
hand, for asymptotically large S and fixed Lc/l > (i.e., fixed iV > and S -^ oo), the number of tethered sites 
N becomes negligible compared to cluster size S. Thus, the lowest energy mode is arbitrarily close to stationarity 
and should behave essentially in the same way as the second lowest energy mode (lowest energy, nontrivial mode) of 
the nontethered case, i.e., as ei ~ g-d^/df (^^^ jg |-]-^g ^q//j dimension and df is the fractal dimension of the cluster). 
Thus, the scaling form of ei is expected to be 

ei ^ S-'^-/'^fF{Ljl) = S-'^-/'^fG{N). (6) 



Since ei ^ as A^ — ^ with fixed S, we must have G{z) ^ as z ^ 0. If we assume analyticity of G{z) for 
small z, in the absence of obvious symmetry requirements which might remove the linear term, it seems reasonable 
to conclude G{z) ~ z for small z. 

It is possible to obtain the limit of z — *■ cxd in such a way as to essentially remove the effect of tethering simultaneously. 
To this end, consider joining many mini-clusters with a given individual h in such a way so as to make the overall 
cluster a fractal with fractal dimension df. In this process the overall cluster length scale Lc increases at a greater 
rate (Lc ~ S^^'^f) than the overall I (/ ~ (S/NY^'^f) (since iV is also growing). Thus, as more mini-clusters are joined, 
the variable Lc/l tends to oo while the tethered hull (S'^f''^^) becomes a negligible part of the entire boundary which 
goes as S (since dh < df). In this case with the diminished importance of tethered hull, the behavior of ei should 
converge toward the lowest nontrivial mode of the nontethered limit once again (ei ~ S~'^"'^'^f). Thus, in this limit 
of z = iV — > oo, the lowest eigenmode is confined entirely to the interior of the cluster, and the spatial extent and the 
energy of the mode is unaffected by further tethering of hull sites. This is equivalent to herding |9|] which is a term 
we had used to describe the saturation effect of the confinement of an eigenmode. Thus, we must have G{z) — > const. 
as A^ — > oo. 

Next we consider all boundary tethering. The no-tethering limit is of course the same as for the hull tethering 
case. As one tethers sites of all boundaries, the lowest-energy mode gets localized in blobs which are unaffected by 
tethering. Moreover, if the ratio Lb/l is fixed at a nonzero value and S — > oo, then each blob behaves like a cluster 
of its own with hull tethering. The expected scaling variable is then the number of tethered sites per blob NL^^' / S 
(which is equal to {Li,/l)'^f as expected). We thus propose a scaling form 

ei - S-'^-Z'^fFiLb/l) = S-'^-/'^'G{N/S''), (7) 

where x = 1 — ydf . 

For z ^ 0, as for the hull tethering, we have G{z) -^ 0. However, the limit of z -^ oo can be achieved by increasing 
the fraction of tethered boundary sites as the cluster size increases. The limit should then correspond to the so-called 
ideal chain flo] with all the boundary sites tethered. Thus, in this limit, N ^ S and ei ^ (InS')^^'''" |jl9] (where 
do is the exponent describing the stretched exponential behavior of the density of states for ideal chains introduced 
in pql). Numerically, do is about 4 in two dimensions and about 6 in three dimensions for the critical percolation 
cluster. Thus, the scaling function G{z) ^ 2;'^™/['^/(i-^)](ln2)~2/''o as z — * oo. 

It is interesting to consider the connection of this behavior to the confinement of the lowest energy mode to the 
small, compact regions of the cluster in this case. These regions are not affected by tethering because of the absence 
of internal boundaries (holes). Their fractal dimension must be essentially equal to the Euclidean lattice dimension 
d. Thus, we denote their average size by Sc, we may expect ei ^ Sc~ where 2 is the walk dimension on a compact 
cluster. In view of the relation ei ^ {\nS)~^^'^° above, this would lead to Sc ~ {\nSY^'^° ~ (In 5)^/^ in both two and 
three dimensions. This type of slow growth of Sc is consistent with our direct observations. 

Fig. 2 and Fig. 3 show the numerical results for the above discussed scaling behavior for the square lattice in two 
dimensions. Fig. 2 is for the hull tethering case while Fig. 3 is for all boundary tethering. Besides the scaling variables 
on the X axis being different, the scaling functions exhibit dramatically different behavior. Though, in both cases, ei 
(and thus the scaling functions) approach zero as z — > 0, the low z behavior is linear for hull tethering while it seems 
to be faster than linear for all boundary tethering. Also while the hull tethering scaling function tends to a finite limit 
as 2 ^ oo (see above discussion), that for all boundary tethering appears to grow unbounded as it should if a; < 1. 
Corresponding results for all boundary tethering for the simple cubic lattice are shown in Fig. 4. Here the reasonable 
data collapsing is achieved for x = 0.4. 

We find, in particular, that hull tethering cannot change the behavior of the lowest energy mode, which remains 
the same as that of the untethered case (ei ^ S~'^"'^'''f). This is because a fully tethered hull confines the mode in 
the interior of the critical percolation cluster which, asymptotically for large clusters, scales with the same fractal 
dimension as the entire cluster itself. Thus, the lowest energy mode results in the vibration of the interior of the 
cluster independent of the external boundary. This inability to change the behavior of the lowest energy mode is 
analogous to the case of an Euclidean cluster as well as to the fractal drums (fractal external boundary and Euclidean 
interior) g). 

In contrast, the tethering of all boundaries leads to confinement of the lowest energy vibrational mode in regions 
of the cluster which are essentially compact, altering its behavior from a power-law dependence on cluster size 
£1 ^ g-duj/df ^Q g^ slow logarithmic dependence on ei '^ {h\S)~'^^'^° . Thus all boundary tethering conspires with the 
complex geometry of the critical percolation cluster (with pockets of compact regions, etc.) to dramatically alter the 
behavior of lowest energy vibrational mode. 



III. SCALING OF THE MAXIMUM IN DENSITY OF STATES 

Let us now consider the vibrational density of states p(e) of the matrix W. For nontethered hmit, p(e) generaUy 
has a power law increase toward the lowest energy (stationary) mode Q. On the other hand, if all boundary sites are 
tethered, we recover the ideal chain result of a stretched exponential decrease toward the lowest energy mode ei [ p^ . 
Thus, when progressively more of all boundary sites are tethered, a crossover between the two limits occurs. For an 
intermediate fraction / of tethered sites, we generally observe a maximum in the density of states as shown in Fig. 5(a) 
for the square lattice. Starting from larger e and moving toward e = 0, p increases initially because the nontethered 
modes for the range of e are too localized to be affected when tethers are added, while for much smaller e the tethers 
begin to drag down the density of states. (Remember that tethers correspond to traps and low e corresponds to the 
long time survival of a diffusing particle.) Thus a peak in p(e) occurs, say, at gp, which gradually increases from for 
/ = {ant limit) toward f — I (ideal chain limit) as the tethering fraction / increases. 

The existence of the maximum in the density of states has implications on the possible resonant behavior in response 
to the external source of vibrational energy. See Leibig H for a discussion of this point for the case of weakly disordered 
network of Hookian springs. 

For the hull tethering case, the situation is similar for finite S. In the asymptotic large S limit, however, hull 
tethering cannot affect the density of states in qualitative manner no matter how large the tethering fraction / may 
be because the hull becomes increasingly a negligible fraction of all boundaries. This implies that, for large S, the 
region near e = must always witness power-law increasing p{e). However, the finite S effect often masks this 
asymptotic behavior so that for most sizes of clusters numerically generated and for most tether fraction /, we do 
observe a maximum in the density of states similar to the all boundary tethering case (see Fig. 5(b)). 

The crossover of £p as the tether fraction / is varied is then expected to obey the same type of scaling as for the 
lowest energy mode ei. For hull tethering, 



while, for all boundary tethering. 



S-d^-/dfGp{N), (8) 



S-d^^/'^fGpiN/S^), (9) 



with the same value of x as in Eq.(|7[). 

These scaling laws have been tested on the same critical percolation clusters as for the lowest energy mode scaling. 
The numerical results, shown in Fig. 6 for the hull case and Fig. 7 for all boundary tethering, are in good agreement 
with the scaling laws above. The corresponding results for the simple cubic lattice in three dimensions are shown in 
Fig. 8, again with reasonably good agreement when a choice oi x = 0.4 is used. 

IV. SPATIAL EXTENT OF THE LOWEST ENERGY MODE 

The difference between hull and all boundary tethering is dramatic also in the spatial characters of the normal 
modes. We show in Fig. 9 and Fig. 10 the amplitude maps of typical lowest energy modes (corresponding to ei) for 
hull tethering and all boundary tethering, respectively. For hull tethering, starting from low tethering fraction / (or 
small N) and increasing / initially reduces the amplitude of vibration of the sites in the vicinity of the tethered hull 
sites. This causes a decrease in the wavelength and consequently an increase in ei (also of the scaling function F{z)). 
However after a certain number of hull sites are tethered, the amplitudes of the mode at all sites near the hull become 
attenuated and the region of large amplitudes becomes well confined to the interior of the cluster. Beyond this point, 
increasing N does not further affect the spatial structure of the normal mode as the sites with appreciable amplitudes 
are already deeply in the interior of the cluster. This saturation behavior may be likened to herding of livestock into 
a safe, fenced haven, and thus we may call it herding of normal modes. Herding is reflected in the saturation behavior 
of F{z) for z tending to infinity. In the examples shown in Fig. 9, the participation rations indicate that only about 
1% of the sites contribute substantially to the mode for all shown values of the tethering fractions. 

For all boundary tethering, increasing / serves to progressively confine the (originally most extended, untethered) 
normal mode to sites which are interior to both the external and internal boundaries, i.e., in the blobs of the percolation 
cluster. Since typically most blobs are very small, this results in a much sharper decrease in the spatial extent of 
the normal mode, as seen in Fig. 10. Even for relatively small /, the effect of confinement can be nearly complete in 
this case, and further tethering may force the mode to jump around to distant locations which may allow the largest 
spatial extent. This behavior may be likened to chasing wild animals in a hunting expedition, and thus we may call 
it hunting of normal modes. It is interesting to note that, while the largest modes without tethering are those of 



small e, they are also the ones most affected by tethering (particularly all boundary tethering), and thus after some 
tethers have been put in place, these modes are not necessarily the spatially largest ones any longer. Indeed those 
modes shown in Fig. 10 have main contributions only from less than 0.1% of the sites according to the values of the 
participation ratios. 

The above observations on the spatial structure of normal modes under boundary tethering may have technological 
implications. For example, a soft glassy material may need to be clamped only at relatively few locations on its 
external boundary to fully confine its lowest energy normal mode to its interior via herding. 

V. CONCLUSION 

The present work reveals a number of interesting features of the elastic properties of a fractal with tethered 
boundaries. The vibrational modes and consequently the vibrational density of states for both hull and all boundary 
tethering are dictated by regions of clusters which are tether free. Despite the structural similarity between the hull 
and internal boundaries, the effect of tethering depends greatly on whether only the hull is tethered or all boundaries 
are tethered. Tethering of hull sites confines the low-energy vibrational modes in the interior of the cluster and since 
the interior of the cluster is itself a fractal with the fractal dimension same as that of the entire cluster, the leading 
behavior of the vibrational density of states for fully tethered hull must asymptotically be the same as that of the 
untethered cluster. On the other hand, for full tethering of all boundaries, the low-energy modes are confined in 
the compact regions of the cluster which, because of their small spatial extent, give rise to faster than exponential 
decrease in the density of states in the same limit. 

Keeping in mind the fact that Euclidean clamped boundaries do not have an appreciable effect on the density 
of states, our results demonstrate the potential of tethering fractal boundaries to attenuate harmonic excitations. 
Also, the present problem is equivalent to that of diffusion in the presence of permanent traps, a fact that was 
taken advantage of for the numerical part of this work. Since the latter serves as a model for diffusion controlled 
reaction/absorption in the presence of immobile reactants/absorbents, there may also be applications in the areas 
of drug reaction/absorption, etc. Finally, because of the similar mathematical formulation of quantum mechanical 
localization and hopping transport problems Q , we expect that the same techniques will be useful in studying those 
problems and that some of the current results may have direct analogs in them. 
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FIG. 1. The average distance I between tethered sites, average diameter of the cluster Lc, and the average blob diameter Li,, 
for (a) all boundary tethering and (b) hull tethering. 

FIG. 2. Scaling of the lowest energy mode for hull tethering on square lattice. Scaling variable reduces to A'^, which is the 
total number of tethered hull sites. 

FIG. 3. Scaling of the lowest energy mode for all boundary tethering on square lattice. The scaling variable is N/ S'^ where 
X is numerically about 0.54. 

FIG. 4. Scaling of the lowest energy mode for all boundary tethering on simple cubic lattice. The scaling variable is N/S^ 
where x is numerically about 0.4. 

FIG. 5. (a) Density of energy eigenvalues are shown for all boundary tethering of 5000-site clusters at the tethering fraction 
of 0.1, 0.4, and 0.6. (b) shows the same for hull only tethering at fraction 0.3, 0.5, and 0.9. 

FIG. 6. Scaling of the peak in the density of states for hull tethering on square lattice. 

FIG. 7. Scaling of the peak in the density of states for all boundary tethering on square lattice. 

FIG. 8. Scaling of the peak in the density of states for all boundary tethering on simple cubic lattice. 

FIG. 9. Spatial extent of the lowest energy mode is illustrated by plotting the amplitudes of the mode for a 50000 site cluster 
on the square lattice for progressively increasing hull tethering fraction. Grey scale is used where lightest grey is used for the 
lowest (absolute value of) amplitude and black for the highest, in logarithmic scale, (a) corresponds to random tethering of 
20% (upper left), (b) to 40% (upper right), (c) to 70% (lower left), and to 100% (lower right). Note that with zero tethering, 
the amplitude will be completely uniform. 

FIG. 10. Same as Fig. 9 but for the all boundary tethering case. 
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